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Table 1. ALFA ZOA Deep Survey: First Results Catalog. 


ALFAZOA 

RA 

Dec 

1 

b 

Fhi 

Vhel 

bFso 

W 20 

Dlg 

IoqMhi 


(J2000.0) 

(J2000.0) 

n 

n 

(Jy km ) 

(km ) 

(km ) 

(km ) 

(Mpc) 

(MJ 

J0555 + 1531 

05 55 43 : 

+15 31 54 : 

192.73 : 

-4.84 : 


6620 : 





J0556 + 1520" 

05 56 11 

+15 20 41 

192.95 

-4.84 

2.5 + 0.2 

7870 + 1 

210+3 

222+4 

116 

9.9 

J0556+1450 

05 56 30 

+14 50 34 

193.43 

-5.02 

4.7+0.3 

6774+6 

326+12 

367+19 

100 

10.0 

J0556+1513 

05 56 54 

+15 13 22 

193.14 

-4.74 

0.6 + 0.1 

7999 + 6 

137+13 

184+19 

118 

9.3 

J0556+1448 

05 56 55 

+14 48 25 

193.51 

-4.95 

0.6 + 0.1 

5474+4 

152+8 

170+11 

81 

8.9 

J0557+1450 

05 57 02 

+14 50 16 

193.49 

-4.91 

0.8 + 0.1 

8070 + 7 

176+14 

230+21 

119 

9.4 

J0558 + 1708" 

05 58 08 

+17 08 51 

191.61 

-3.53 

1.0 + 0.1 

4966+4 

164+8 

193+12 

73 

9.1 

40558+1557'= 

05 58 48 

+15 57 31 

192.73 

-3.98 

1.1 + 0.1 

6966 + 3 

168+5 

190 + 8 

103 

9.4 

J0558+1656 

05 58 52 

+16 56 03 

191.89 

-3.48 

0.8 + 0.1 

8222 + 2 

114+5 

141 + 7 

121 

9.5 

J0559+1636 

05 59 29 

+16 36 40 

192.24 

-3.52 

1.7+0.1 

8039 + 2 

327+5 

352 + 7 

118 

9.8 

40559+1543'= 

05 59 31 

+15 43 43 

193.02 

-3.95 

2.1 + 0.2 

6529 + 3 

74+6 

112 + 9 

96 

9.7 

40559+1532'= 

05 59 49 

+15 32 23 

193.22 

-3.98 

6.6 + 0.3 

5451 + 1 

188+2 

203 + 3 

80 

10.0 

40600+1644'= 

06 00 31 

+16 44 39 

192.25 

-3.24 

0.6 + 0.1 

8109 + 3 

26+6 

102 + 9 

119 

9.3 

40600+1654 

06 00 37 

+16 54 34 

192.12 

-3.13 

1.5 + 0.1 

7006 + 2 

251+4 

265 + 6 

103 

9.6 
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RA 

(J2000.0) 

Dec 

(J2000.0) 

1 

n 

b 

n 
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(Jy km ) 
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(km ) 

W"50 

(km ) 

W 20 

(km ) 

Dlg 

(Mpc) 

logMui 

(MJ 

J0600 + 1642" 

06 

00 

39 

+ 16 42 40 

192.31 

-3.21 

l.O+O.l 

8170 + 2 

139+4 

157+6 

120 

9.5 

J0600+1640^ 

06 

00 

48 

+ 16 40 41 

192.32 

-3.21 

0.9+0.1 

8206+4 

72 + 8 

148+11 

121 

9.5 

J0600+165P 

06 

00 

52 

+16 51 08 

192.20 

-3.11 

4.5+0.4 

6907+5 

211 + 10 

284+15 

102 

10.0 

J0600+1529 

06 

00 

59 

+15 29 16 

193.40 

-3.76 

1.2+0.1 

5579 + 3 

165 + 5 

194+8 

82 

9.3 

J0601 + 1540 

06 

01 

02 

+ 15 40 14 

193.25 

-3.66 

1.1+0.1 

5509 + 1 

214+3 

220+4 

81 

9.2 

40601 + 1709^ 

06 

01 

11 

+ 17 09 18 

191.97 

-2.89 

2.5+0.2 

6863 + 2 

283+4 

297+6 

101 

9.8 

40601 + 1620^ 

06 

01 

34 

+16 20 20 

192.73 

-3.22 

2.9+0.3 

4010 + 2 

132 + 5 

159+7 

59 

9.4 

40602 + 1558 

06 

02 

03 

+15 58 14 

193.11 

-3.30 

1.1+0.1 

5397+3 

135 + 7 

163+10 

80 

9.2 

40602 + 1450^ 

06 

02 

26 

+ 14 50 51 

194.13 

-3.77 

149 + 0.9 

769+0 

233+0 

250+0 

11 

9.7 

40602 + 1449 

06 

02 

27 

+ 14 49 44 

194.15 

-3.77 

0.3+0.1 

18668 + 5 

61 + 10 

90 + 16 

275 

9.7 

40602 + 1438 

06 

02 

37 

+ 14 38 19 

194.34 

-3.83 

0.5+0.1 

14213 + 5 

152 + 10 

179+15 

209 

9.7 

40603+1601 

06 

03 

05 

+16 01 08 

193.19 

-3.06 

1.9+0.1 

6332 + 2 

233+5 

261 + 7 

93 

9.6 

40603+1607^ 

06 

03 

06 

+16 07 29 

193.09 

-3.00 

3.6+0.3 

5444+3 

310+6 

337+10 

80 

9.7 

40603+1642 

06 

03 

14 

+ 16 42 02 

192.61 

-2.69 

0.7+0.1 

8295+4 

73 + 8 

118+12 

122 

9.5 
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ALFAZOA 

RA 

(J2000.0) 

Dec 

(J2000.0) 

1 

n 

h 

n 

Fhi 

(Jy km ) 

Vhel 

(km ) 

W"50 

(km ) 

W 20 

(km ) 

Dig 

(Mpc) 

logMni 

(Me) 

J0605 + 1526" 

06 05 54 

+15 26 07 

194.03 

-2.75 

3.9+0.3 

5452 + 2 

153+4 

180+6 

80 

9.8 

J0606+1459 

06 06 47 

+ 14 59 08 

194.52 

-2.78 

1.3+0.1 

3911 + 2 

146+5 

162 + 7 

58 

9.0 

J0607+1505 

06 07 17 

+ 15 05 55 

194.48 

-2.62 

1.9+0.1 

5376 + 3 

126+7 

243+10 

79 

9.4 

J0607+1606 

06 07 29 

+ 16 06 00 

193.63 

-2.09 

2.6+0.1 

5637+3 

130+5 

238+8 

83 

9.6 

J0607+1615 

06 07 33 

+ 16 15 00 

193.51 

-2.01 

1.2+0.1 

5515 + 6 

153+11 

268+17 

81 

9.3 

J0607+1543 

06 07 39 

+ 15 43 45 

193.97 

-2.24 

0.4+0.1 

12862 + 5 

80 + 9 

105 + 14 

189 

9.6 

J0607+1507 

06 07 43 

+15 07 43 

194.51 

-2.52 

1.1+0.1 

5436 + 5 

170+10 

216+14 

80 

9.2 

40607+1607^ 

06 07 45 

+ 16 07 00 

193.65 

-2.03 

4.4+0.1 

5573 + 2 

188+4 

253+5 

82 

9.9 

J0607+1611a 

06 07 46 

+ 16 11 00 

193.59 

-1.99 

3.3+0.1 

5539 + 2 

127+4 

203+5 

82 

9.7 

J0607+16116 

06 07 46 

+ 16 11 00 

193.59 

-1.99 

0.7+0.1 

5217+2 

187+5 

196+7 

77 

9.0 

J0608+1527 

06 08 50 

+ 15 27 27 

194.35 

-2.12 

1.0+0.2 

13787+5 

150+10 

182 + 15 

203 

10.0 

J0610+1709 

06 10 20 

+ 17 09 11 

193.04 

-0.99 

2.4+0.2 

13453+3 

215+6 

249+9 

198 

10.3 

J0610+1704 

06 10 57 

+ 17 04 55 

193.17 

-0.89 

0.5+0.1 

13558+4 

143+8 

161 + 11 

200 

9.7 

40611 + 1552^ 

06 11 07 

+15 52 16 

194.25 

-1.44 

12 + 0.5 

5443 + 3 

387+7 

484+10 

80 

10.3 
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ALFAZOA 

RA 

(J2000.0) 
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(J2000.0) 
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n 
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(Jy km ) 
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(km ) 

W"50 
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(MJ 

J0611 + 1553 

06 

11 26 

+15 53 00 

194.28 

-1.37 

0.9+0.1 

5589 + 3 

108+6 

141+9 

82 

9.1 

J0611 + 1643" 

06 

11 31 

+ 16 43 18 

193.55 

-0.95 

0.9+0.1 

5241 + 2 

89+4 

101 + 5 

77 

9.1 

J0612 + 1636 

06 

12 37 

+16 36 29 

193.78 

-0.77 

0.6+0.1 

5146 + 3 

184+5 

196+8 

76 

8.9 

J0613+1545 

06 

13 41 

+15 45 07 

194.65 

-0.96 

1.4+0.2 

13307+3 

202+6 

215+9 

196 

10.1 

J0613+1707 

06 

13 59 

+17 07 41 

193.48 

-0.23 

0.5+0.1 

16308 + 7 

64+14 

148+20 

240 

9.8 

J0614+1641 

06 

14 23 

+ 16 41 30 

193.91 

-0.36 

0.6+0.1 

13366+4 

256+8 

267+12 

197 

9.8 

J0615 + 1547 

06 

15 08 

+ 15 47 38 

194.78 

-0.63 

0.3+0.1 

13062+4 

47+8 

64+12 

192 

9.3 

J0615 + 1705 

06 

15 40 

+ 17 05 47 

193.70 

0.10 

0.5+0.1 

4841 + 5 

91 + 10 

139+14 

71 

8.8 

J0616+1639 

06 

16 00 

+ 16 39 58 

194.12 

-0.03 

0.3+0.1 

2144+4 

38 + 7 

65 + 11 

32 

7.9 

J0616+1526 

06 

16 03 

+15 26 10 

195.20 

-0.61 

0.7+0.1 

11735 + 5 

178 + 10 

213+16 

173 

9.7 

J0616+164P 

06 

16 10 

+ 16 41 54 

194.11 

0.02 

6.8+0.3 

5555 + 3 

210+5 

272+8 

82 

10.0 

J0616+1704 

06 

16 25 

+ 17 04 56 

193.80 

0.25 

0.8+0.1 

5040 + 2 

170+4 

190+6 

74 

9.0 

J0616+1705 

06 

16 34 

+ 17 05 12 

193.81 

0.29 

1.3+0.1 

5704+1 

56 + 2 

86 + 384 

9.3 

9.3 

J0616+1626 

06 

16 45 

+16 26 54 

194.39 

0.02 

1.2+0.1 

11584+3 

386+5 

397+8 

170 

9.9 






ALFAZOA 

RA 

Dec 

1 

h 


(J2000.0) 

(J2000.0) 

n 

n 

J0616 + 1536 

06 16 58 

+15 36 41 

195.15 

-0.33 

J0617+1655 

06 17 15 

+16 55 43 

194.03 

0.36 

40617+1648'= 

06 17 19 

+16 48 26 

194.14 

0.31 

J0617+1644 

06 17 43 

+16 44 55 

194.24 

0.37 


J0618 + 1430 06 18 18: +14 30 53: 196.27: -0.57: 


ontinued 


Fhi 

(Jy km ) 

Vhel 

(km ) 

VFso 

(km ) 

4+20 

(km ) 

Dlg 

(Mpc) 

IoqMhi 

(Me) 

1.3 + 0.1 

11761 + 2 

148+4 

159 + 5 

173 

10.0 

1.2 + 0.1 

5373+4 

369 + 7 

389 + 11 

79 

9.3 

1.4+0.2 

5696+3 

111 + 5 

123 + 8 

84 

9.4 

0.3 + 0 

5313+3 

22+6 

74+9 

78 

8.7 


13667 : 
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ABSTRACT 


The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Deep 
Survey is the deepest and most sensitive blind Hi survey undertaken in the 
ZOA. ALFA ZOA Deep will cover about 300 square degrees of sky behind 
the Galactic plane in both the inner (30° < I < 75°; 6 < |2°|) and outer 
(175° < I < 207°; —2° < b < +1°) Galaxy, using the Arecibo Radio Tele¬ 
scope. First results from the survey have found 61 galaxies within a 15 square 
degree area centered on I = 192° and b = —2°. The survey reached its expected 
sensitivity of rms = 1 mJy at 9 km s“^ channel resolution, and is shown to be 
complete above integrated flux, Fhi = 0.5 Jy km s“F The positional accuracy 
of the survey is 28" and detections are found out to a recessional velocity of 
nearly 19,000 km s“^. The survey conhrms the extent of the Orion and Abell 
539 clusters behind the plane of the Milky Way and discovers expansive voids, 
at 10,000 km s“^ and 18,000 km s“^. 26 detections (43%) have a counterpart in 
the literature, but only two of these have known redshift. Gounterparts are 20% 
less common beyond Vhei = 10,000 km s“^ and 33% less common at extinctions 
higher than Ab = 3.5 mag. ALFA ZOA Deep survey is able to probe large scale 
structure beyond redshifts that even the most modern wide-angle surveys have 
been able to detect in the Zone of Avoidance at any wavelength. 

Subject headings: catalogs — galaxies: distances and redshifts — large-scale structure 
of universe — methods: data analysis — methods: observational — surveys 
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1. Introduction 

It has been suggested (Loeb & Narayan 2008) that undiscovered mass behind the Milky 
Way may explain the discrepancy between the cosmic microwave background dipole and the 
mass density dipole (Erdogdu et al. 2006), which drives the peculiar velocity of the Local 
Group. Neutral hydrogen (H l) surveys of the Zone of Avoidance (ZOA) can successfully 
map large scale structure in the Universe behind the Milky Way (ALFA ZOA; Henning et 
al. 2010, HIZOA; Henning et al. 2005). While shallow Hi surveys have uncovered a large, 
spiral galaxy at ~ 3 Mpc (Kraan-Korteweg et al. 1994), and more sensitive surveys have 
discovered nearby dwarf galaxies (McIntyre et al. 2011, Massey et al. 2003, Begum et al. 
2005), Hi ZOA surveys have shown that there are no hidden Local Group galaxies with 
neutral hydrogen mass, Mhi, greater than Mhi ~ 10® Mq in the southern sky (Henning et 
al. 2000) and Mhi ~ 10^ Mq in the northern sky (Henning et al. 1998). A nearby massive 
spiral (IO^^Mq) behind the Milky Way can now largely be ruled out as a way to recover the 
mass dipole vector, though H l surveys are not sensitive to giant elliptical galaxies (Roberts 
& Haynes 1994). A galaxy cluster (IO^^Mq) at 20 Mpc could affect the mass dipole vector, 
and results from HIZOA in the southern sky are forthcoming. Most of the northern sky is 
only starting to be surveyed deeply enough to detect a cluster at 20 Mpc (EBHIS; Kerp et 
al. 2011). 

More recent results have shown that our knowledge of the bulk flow of peculiar 
velocities within 60 h“^ Mpc (h~^ = Hq/IOO, Hq is the Hubble constant in km s“^ Mpc“^) 
is complete enough to be constrained by cosmic variance and so no new galaxy discoveries 
within this distance will result in a significant correction to the dipole vector (Watkins et 
al. 2009). Asymmetry in large scale structure beyond 100 h”^ Mpc remains as a likely 
explanation (Lavaux et al. 2010), though even the most modern all sky redshift surveys 
retain a Zone of Avoidance at these distances (Huchra et al. 2012), forcing bulk motion 
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studies to “fill in” the ZOA by cloning galaxies above and below the Plane (Lavaux & 
Hudson 2011). 

The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Deep Survey is a 
blind survey for galaxies in the ZOA out to distances of 200 h“^ Mpc. The full survey will 
cover nearly 300 square degrees behind both the inner and outer Galaxy regions visible 
with the Arecibo Telescope, with a sensitivity of 1 mJy per beam (at 9 km s“^ velocity 
resolution). The sensitivity and depth of ALFA ZOA will probe large scale structure out to 
200 h“^ Mpc in an area of sky that intersects with known overdensities: the Orion Cluster, 
Abell 539, Abell 634, and GIZA J0603.8+2939 as well as Cl, C5, CIO, C21, and C31 
(Erdogdu et ah 2006). It will also probe the Orion, Canis Major, and Lepus Voids. The 
survey will detect Hi to greater depth than recent, large Hi surveys (ALFALFA; Haynes 
et ah 2011, HIPASS; Zwaan et al. 2005), allowing ALFA ZOA to probe environments 
at farther distances. The low mass end of the Hi mass function (HIMF) has been shown 
by some studies to have a signihcantly steeper slope when measured at high sensitivity at 
distances outside the local Universe (Arecibo Galaxy Environment Survey “AGES”; Davies 
et al. 2011, Arecibo Ultra Deep Survey “AUDS”; Freudling et al. 2011). If conhrmed, 
this could have a major impact on the amount of hydrogen gas available for ongoing star 
formation and galaxy evolution. ALFA ZOA Deep is similar in sensitivity to AGES and 
surveys 10^ times more area than AUDS, putting it in a unique position to check the latest 
HIMF results as an unbiased, blind survey. 

This paper describes the observation, data reduction, and source detection and 
parameterization techniques for the ALFA ZOA Deep Survey. It presents a catalog of the 
hrst results of the survey, coming from a 15 square degree region, and uses these results to 
analyze survey performance. Section [2] describes the observation technique of the survey. 
Section [3] discusses the data reduction process. Section 0] describes the survey search 


method and source parameterization. Section [5] presents the hrst results catalog. Section 
E] discusses survey accuracy, sensitivity, and effectiveness. Section [7] discusses large scale 
structure uncovered by the survey. Section [H] is the conclusion. 


2. Survey Design 
2.1. Receiver 

ALFA ZOA utilizes the ALFA receiver on the 305m Arecibo Radio Telescope in 
Puerto Rico. The ALFA receiver consists of seven independent beams with two orthogonal 
linear polarizations each, allowing the survey to cover the same area seven times faster 
than a single beam receiver. The center beam is surrounded by six outer beams in a 
hexagonal pattern. The receiver covers frequencies 1-2 GHz and the FWHM at 1.4 GHz is 
approximately 3.4' per beam. The gain of the seven ALFA beams ranges between 8.5 and 
11 K Jy“^, and the system temperatures range between 26 and 30 K. 

Observations are recorded in one-second integrations using the Mock Spectrometer, 
which performs an “on the fly” Fast Fourier Transform each time the voltage is sampled. 
The digital spectrometer covers 300 MHz from 1225 MHz to 1525 MHz using two 
overlapping 172 MHz sub-bands centered on 1300 MHz and 1450 MHz. The sub-bands 
overlap between 1364 MHz and 1386 MHz so that there is no loss of sensitivity due to roll 
off at the edges of the bandpass. Each sub-band is divided into 8192 channels, producing 
a spectral resolution of 21 kHz. This equates to a velocity resolution of ~ 4.5 km s“^ for 
neutral hydrogen emission. 
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2.2. Observations 

The ALFA ZOA Deep survey covers nearly 300 square degrees through both the inner 
(30° < I < 75°; fe < |2°|) and outer (175° < I < 207°; -2° < 6 < +1°) Galaxy. ALFA ZOA 
takes data simultaneously with a survey for pulsars in the Milky Way (PALFA; e.g. Cordes 
et ah 2006) and a Galactic radio recombination line survey (SIGGMA; Liu et al. 2013), 
as well as several SETI groups that receive data from the ALFA receiver (Astropulse; Von 
Korff et al. 2013, SETI@home; Anderson et al. 2002, Serendip V.v; Gobb et al. 2000). In 
particular, the hrst fast radio burst detected with an instrument other than the 13-beam 
receiver of the Parkes Radio Telescope (Spitler et al. 2014) was discovered by PALFA using 
the same data taken by ALFA ZOA for this paper. 

Observations in the outer Galaxy are controlled by ALFA ZOA, and the inner Galaxy 
is controlled by PALFA. The setup of observations is slightly different between the inner 
and outer regions and this is discussed in the next section. The hrst results presented in 
this paper come from a completely surveyed, 15 square degree area in the outer Galaxy 
centered on 1 = 192° and b = -2°. The area was chosen to intersect large scale structure 
predicted from known structure above and below the Galactic Plane. The dimensions are 
330' X 164', ranging across right ascensions 05:55:30 to 06:18:20 and declinations -1-14:30:00 
to -1-17:12:00. Observations were conducted from December 2010 to March 2012. Figure 
[T] shows a map of the survey area in the outer Galaxy including pointings that have been 
observed so far. 


2.3. Tiling Pattern 

Gommensal observations are an incredibly efficient use of telescope time (e.g. hve other 
projects listed above take data simultaneously with ALFA ZOA) but require compromises 
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Fig. 1.— Sky coverage of the survey in the Outer Galaxy in Galactic coordinates. Each black dot is a 
central beam pointing. The area inside the dashed, red box is the focus of the First Results Gatalog. The 
area inside the solid, green box is the focus for the full survey. 


between commensal partners. While drift scans have been shown to produce superior 
baselines for spectral line surveys (Briggs et al. 1997), the ALFA ZOA Deep Survey is 
constrained to tracking single pointings due to commensal obligations. One method to 
bandpass correct a single pointing is self-subtraction of a bandpass by its median hlter, 
though this was shown to be inadequate for use as a data reduction procedure in a blind H l 
galaxy survey (McIntyre 2013a). As such, the survey uses the position switching method 
for bandpass correction, a non-standard observing mode for extragalactic ALFA surveys, 
requiring ALFA ZOA Deep to observe position switched pairs and develop unique software 
for data reduction. Position switching uses a total power on-off technique to remove the 
bandpass from a source (ON) by subtracting the bandpass from a position off-source (OFF). 
ALFA ZOA uses the geometry of the ALFA pattern to generate a list of pointings that 
most efficiently covers the survey area and is designed so that several different pointings 
can be tracked over the same path of alt-az coordinates in order to produce flat baselines. 

The survey tiling pattern is made from clusters of three beam patterns, the centers of 
which are separated by one FWHM along an axis tilted 19° from north-south, as shown 
in Figure [21 Clusters fully cover their area and £t into surrounding clusters in order to £11 
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Fig. 2.— Tiling pattern used by the ALFA ZOA Deep Survey shown for a small region of sky. Clusters 
of three pointings fully cover their area and fit into surrounding clusters. The seven beams of the ALFA 
pattern are included for two of the clusters to demonstrate how the geometry of the ALFA beams is used to 
cover the survey area. One cluster shows the beams coded separately for the three different pointings. 

in the survey area without Nyquist sampling. The axis of the ALFA receiver is rotated 
counter-clockwise, 19° from north-south, so that adjacent clusters to the east and west are 
centered on the same declination roughly 60 seconds of right ascension apart. This allows 
the survey to integrate one pointing for 270 seconds, for instance, and have time to slew 
to another pointing 300 seconds of right ascension away so that it can then track over the 
same alt-az path (i.e. create an ON-OFF pair). Every two observations must be observed 
in pairs in this way in order to apply bandpass correction, and the full survey area is hlled 
out by leapfrogging through this tiling pattern. 

The observing setup is slightly different between the inner and outer Galaxy regions 
because of commensal obligations. In the inner Galaxy, each pointing is integrated for 270 
seconds on sky and then 6 seconds on the ALFA high flux noise diode for calibration. A 
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single ALFA beam pattern covers a 60 square arcminute area, and about 10 pointings can 
be observed in an hour. At this rate, the survey covers approximately 10 square degrees per 
60 hours of telescope time, though the area is not covered at the Nyquist sampling rate. 

Outer Galaxy observations are taken for 180 second integrations and 6 seconds on 
the high flux noise diode. Each pointing has two Nyquist counterparts, each 1/2 FWHM 
away on either side, perpendicular to the axis of the three pointing cluster. The effective 
integration time for each spatial position is the same as in the inner Galaxy. However, the 
survey coverage rate is slower by a factor of two because three 180 second integrations are 
taken in the outer Galaxy for every one 270 second integration taken in the inner Galaxy. 
Surveying an area that is Nyquist sampled allows the integrated flux of a detection to be 
fully recovered. 


3. Data Reduction 

3.1. Process 

The ALFA ZOA Deep Survey developed unique software for data reduction in the 
IDL programming language because of the total power on-off observing technique. Each 
one-second integration of a telescope pointing needs to be identihed with its off-source 
counterpart, within a small angular separation, 66. In an effort to dehne an upper limit 
for 66, several observations were taken during June 2010 to test the effect of ON-OFF 
alt-az positional difference on the quality of spectral baselines. The results from this test 
show that the upper limit is 66iim = 1.7' (McIntyre 2013b), or 1/2 FWHM. If a one-second 
integration does not have an off-source integration within 66iim, then it is rejected. For an 
ON that has at least one OFF within 66iim, then the OFF with the smallest 66 for each 
ON integration is used to create a reduced bandpass, (ON-OFF)/OFF. This is divided 
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by the median of the OFF bandpass for normalization. All the reduced spectra from one 
pointing (maximum 270 in the inner Galaxy, 180 in the outer) are median-combined to 
reject spurious RFI and recover the total integration time of the pointing. Most ON-OFF 
pairs have integrations rejected due to separation larger than SOum, though it is very rare to 
have more than 20 integrations rejected. The noise diode from the OFF is used to calibrate 
from receiver units into temperature, and the gain, G(beam, ZA, polarization), converts 
temperature into janskys. The frequency is resampled with a linear interpolation from 
topocentric to barycentric and converted relativistically into heliocentric velocity using 
the optical velocity convention. The two polarizations are averaged and the spectrum is 
Hanning smoothed, which mitigates ringing effects and lowers noise by a factor of \/2, but 
worsens spectral resolution by a factor of two. 

The hnal, reduced spectrum for a 270 second integration has a noise level of 1 mJy 
per channel at a velocity resolution of 9 km s“^. The velocity range of the high frequency 
sub-band is cut to heliocentric velocities, Vhei = -2,000 km s“^ to 12,000 km s“^ and a 
third order polynomial baseline is auto-htted and subtracted over this spectrum. The low 
frequency sub-band is cut to a velocity range of Vhei = 10,000 km s“^ to 21,000 km s“^ 
and also fit by a third order polynomial. The bandpass at velocities beyond 21,000 km s“^ 
contains a signihcant amount of RFI and has not been incorporated into the survey at this 
point. The reduced spectra are then made into data cubes of right ascension, declination, 
and velocity by the program Gridzilla (Barnes et al. 2001). 

Gridzilla is part of the AIPS-I--I- software package, and is used to create a data cube by 
spatially gridding spectra. For ALFA ZOA Deep, Gridzilla is set to clip spectra between 
-50 mJy and 500 mJy, and make a cube containing 1' x V pixels. Each input spectrum is 
assigned to all pixels within a 4' diameter of the spectrum’s coordinates using a Top-Hat 
kernel. In pixels where multiple spectra contribute they are median-combined, weighted by 
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the angular distance of the pixel from each spectrum’s original coordinates using a beam 
power pattern with FWHM = 3.4'. The end result is a 3-dimensional data structure with 
axes of right ascension, declination, and heliocentric velocity. Position-velocity planes can 
be iteratively viewed to search for galaxy detections. 


4. Search Method and Source Parameterization 

The data cube is searched for detections by visually inspecting over the usable velocity 
range using the visualization tool. Karma KVIS (Gooch 1996). This is done by examining 
image planes in position-velocity slices and using a greyscale to represent flux intensity. The 
data cube is searched by three authors (two authors and a UNM graduate student in the 
case of this First Results paper), who look for prohle shapes consistent with known galaxies, 
e.g. Gaussian, double-horn, etc. Lists of galaxy candidates are prepared independently by 
each searcher. These lists are compared with the astronomical coordinate comparison tool, 
Starlink TOPGAT (Tool for OPerations on Gatalogues and Tables) version 3.9 (Taylor 
2005), in order to match sources with RA, Dec, and velocity coordinates within tolerances 
of 5', 5', and 300 km s“^, respectively. The matched lists are adjudicated by a separate 
author who re-examines the position of each candidate source and either accepts it into 
the catalog or rejects it as a false detection. These adjudicated sources are adopted into a 
working catalog and each source is parameterized using MBSPEGT in the software package, 
MIRIAD (Sault et ah 1995). 

For unresolved sources, a spectral profile is created from the weighted emission inside 
a 5' X 5' box surrounding the position of the source. The emission is weighted by spatially 
fitting a Gaussian to a moment 0 map integrated over a user-defined range of velocities, 
and a prohle is created as shown in Figure [3l In an iterative process, the prohle shape is 
visually inspected and a new velocity range is chosen for the prohle (shown as the vertical 
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dotted lines in Fignre[3]). The profile window is also nsed to create a mask for antofitting 
a polynomial to the baseline within 2000 km s“^ of the prohle, nsnally to hrst order bnt 
up to fourth order in rare cases. The integrated flux, Fhi, is computed by subtracting 
the polynomial £t and integrating the flux inside the prohle window. The open circles on 
the outer edges of the prohle in Figure [3] represent the width-maximizer velocity widths, 
IFso and W 20 , at 50% and 20% peak hux, respectively. The closed circle at the top is the 
location of the peak hux. The x near the center of the prohle indicates a width-minimizer 
parameterization, which is not used for ALFA ZOA. The central heliocentric velocity, Vheh 
is taken as the midpoint of the IF 50 value. 


10617+1648 



Fig. 3.— Hi profile of ALFA ZOA J0617+1648, created from MIRIAD output. The vertical dotted lines 
are user defined and create a profile window that source parameters are extracted from. The open circles on 
the outer edges of the detection represent the W 20 and velocity widths, and the closed circle at the top 
is the location of the peak flux. 
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All sources go through the MIRIAD parameterization process as unresolved, but then 
are re-examined to determine the existence of sources with extended emission. A hrst list of 
potentially extended sources is constructed from the size of the Gaussian £t for each source 
as output by MBSPECT. All sources with major axes greater than 3.4' make a candidate 
list of extended sources. Moment 0 maps are then created for each candidate using the 
MIRIAD task, MOMENT. These maps are inspected visually using Karma KPVSLICE 
(Gooch 1996) and a spatial prohle for each source is created by integrating flux across a 
user-defined major axis. These spatial prohles are saved as an array and fed into the IDL 
program GAUSSFIT, which computes a non-linear least-squares Gaussian fit. From this fit, 
the FWHM of the source is determined. 

Every source with FWHM > 3.4' makes it into a shorter candidate list and all of 
these are re-examined again using KPVSLIGE in order to determine an individual aperture 
size to use for each galaxy. KPVSLIGE shows the velocity spectrum from the data cube 
that corresponds with each spatial pixel in the moment map. The user counts how many 
single spatial pixels contain a clear detection in the velocity spectrum, and decides on 
a corresponding aperture box height and width. These new box sizes are then used to 
re-determine source parameterization using MIRIAD, where spectral profiles are created by 
summing flux over all channels and pixels equally inside the aperture and within the profile 
window. The position of the source is determined from the center pixel of the aperture, as 
opposed to the center of a Gaussian £t, but all other source parameterization techniques 
are the same as for unresolved sources. The rms of resolved sources is higher than for 
unresolved sources because of the greater number of pixels summed over. A moment map 
for resolved source, ALFA ZOA J0602-I-1452, is shown in Figure HI The box size of the 
aperture used is 23' x 23' (RA x Dec). 
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Fig. 4. — Moment map of ALFA ZOA J0602+1452. Beam size (FWHM = 3.4') is shown in bottom right 
corner. A clear extension beyond 5' x 5' can be seen on this resolved source. 

5. The ALFA ZOA Deep Survey: First Results Catalog 

The ALFA ZOA Deep Survey: First Results catalog contains 61 galaxies. Profiles are 
shown in Figure O Table 1 presents the following information on the catalog: 

Column (1). ALFA ZOA source name. Sources with extended emission are labeled 
with the ’e’ superscript. 


Column (2). Right Ascension of the htted position in hours, minutes, seconds, J2000 


epoch. 
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Column (3). Declination of the fitted position in degrees, arcminutes, arcseconds, 
J2000 epoch. 

Columns (4) and (5). I and b, Galactic longitude and latitude in degrees, respectively, 
of the fitted position in degrees. 

Column (6). Fhj, integrated flux in Jy km s“^. 

Column (7). Vhei, heliocentric velocity in km s“^. 

Columns (8) and (9). W^q and W 20 , velocity width in km s“^ of the prohle at 50% and 
20 % of the peak flux level, respectively. 

Column (10). distance to the galaxy in Mpc in the Local Group reference frame 

(Courteau & van den Bergh 1999), using Hubble’s Law for cosmological redshift distance 
and taking Hq = 70 km s“^ Mpc“^. 

Column (11). Log Mhi-, logarithm of the total Hi mass in Mq calculated from, 

Mhi = 2.36 X 10^ DIq Ffff, (1) 


where Dig is the distance to the galaxy and F^j is the integrated flux as described above. 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [31 (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
































































































































































Flux Density (Jy) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the prohles is the same 
as in Figured (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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Fig. 5.— Hi profiles of ALFA ZOA detections. The description of the profiles is the same 
as in Figure [3l (Continued) 
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The uncertainties on Fhi, Vheu W 50 , and IT 20 were calculated in the same way as 
Koribalski et ah (2004). The error on the flux integral is 

a{FHi) = 4(^/iV)-i {Sp Fhi 5vY/\ (2) 

where Sp is the peak flux, S/N is the signal-to-noise ratio Sp to cr^Sp), F^i is the integrated 
flux, and Sv is the velocity resolution of the data, 9 km s“^. ct(S'p) is the error in the peak 
flux 

<7{SpY = rms^ + (0.05 SpY- (3) 

The uncertainty in the heliocentric velocity is 


a{vhei)=^{S/N)-^ {P SvY^Y 


( 4 ) 


where 

P = 0.5 (IT 20 - W^5o) (5) 

is a measure of the steepness of the profile edges. The uncertainties in the velocity widths 
are given by 

cr{W 2 Q) = 3 a{vhei), (6) 

^■(fTso) = 2 (T{Vhel)- (7) 

The uncertainties on Dig and Mhi, whose values rely heavily on cosmological assumptions, 
are not calculated. 


Two sources are located on the edge of the cube, J0555+1531 and J0618+1430, and 
so their flux could not be completely recovered. Their positions are listed in the catalog 
followed by a semi-colon in order to indicate that they were not measured through the 
parameterization process described in the previous section. None of their other parameters 
are derived and so are indicated with ellipses. 
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Histograms of parameters from the survey are shown in Figure [ 6 l The distribution of 
heliocentric velocity shows detections out to nearly 19,000 km s“^. The velocity widths 
show detections of dwarf galaxies with H 50 ~ 30 km s“^ as well as large spirals with W^o 
~ 400 km s“^. The distribution of integrated flux ranges from 0.2 Jy km s“^ to 150 Jy km 
s“^. The distribution of mass ranges from Mhi= 10 ^'® to 


5.1. Counterparts 


The NASA/IPAC Extragalactic Database (NED) was searched for extragalactic sources 
within a radius of 2 ' at each galaxy’s position in order to find potential counterparts in the 
literature. NED was also searched for H l extragalactic counterparts within 8 ' and 100 km 
s“^. The counterparts are listed in Table 2, which presents the following information: 

Column (1). ALFA ZOA source name. 

Column ( 2 ). Galactic longitude, I, in degrees. 

Column (3). Galactic latitude, b, in degrees. 

Golumn (4). Foreground extinction, A^, as estimated by Schlafly & Finkbeiner (2011). 

Golumn (5). Primary name of the counterpart as given by NED. 

Golumn ( 6 ). Separation in arcminutes between ALFA ZOA detection and counterpart. 

Golumn (7). Difference in Vhei between ALFA ZOA detection and counterpart in km 


-,-1 


Figure [7] shows the distribution of ALFA ZOA and 2MASS detections within the data 
cube. There are 26 detections (43% of all sources) that have at least one counterpart within 
2 '. Two of these have known redshifts. Some positions have more than one counterpart 
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Fig. 6.— Histograms of HI parameters from the ALFA ZOA Deep Survey. From top to bottom: heliocentric 
velocity, velocity width at 50% peak flux, integrated flux, logarithm of H i mass. 

within 2', and so there are 30 possible counterparts listed in the table. No attempt is made 
to distinguish between multiple counterparts for the same position unless the counterpart 
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Fig. 7.— Distribution of ALFA ZOA (open blue circles) and 2MASS (black dots) detections within the 
survey area. 

has a known Vhei within 100 km of the ALFA ZOA detection. Every Hi source with a 
counterpart has at least one 2 Micron All Sky Survey (2MASS; Skrutskie et ah 2006) galaxy 
associated with it. Figure |8] shows the histogram of ALFA ZOA heliocentric velocities. 
Detections that have a possible counterpart are coded with a diagonal stripe. Counterparts 
from known overdensities can be seen from the Cl cluster at 2000 km s“^, Orion cluster at 
6000 km s“^, Abel 539 at 8000 km s“\ and C21 cluster at 14,000 km s“^ (Erdogdu et ah 
2006). The percent of detections with a counterpart beyond 10,000 km s“^ drops by 20%. 
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Fig. 8.— Histogram of heliocentric velocities for galaxies with no counterpart (solid) and with a counterpart 
(diagonal). 
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5.1.1. Galaxy Classification 


There is no clear bimodal color distribution for galaxies in NIR like there is in optical 
(Jarrett 2000) and so 2MASS colors do not reliably allow for the classihcation of an 
individual galaxy’s morphological type. However, there are parameters other than color 
indices that indicate morphological type, such as mass and size. 
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Fig. 9.— W 50 versus Mhi for ALFA ZOA detections, showing both objects with counterparts (red, open 
circles) and without counterparts (blue, x’s). 

Figure [9] shows velocity width, H 50 , versus Hi mass, Mhi, for ALFA ZOA detections, 
marking objects that have a counterpart and those without a counterpart separately. The 
existence of a counterpart has a strong dependence on velocity width and H l mass. There 
are no counterparts with a velocity width below W^q = 70 km s“^, and every detection with 
a velocity width above HAo = 215 km s“^ has a counterpart. The mean velocity width is 
W 50 = 218 km s“^ for detections with a counterpart and 125 km s“^ for detections with 
no counterparts. The mean Hi mass is log < Mfji > = 9.65 Mq for detections with a 
counterpart versus log < Mhi > = 9.31 Mq for detections without counterparts, over two 
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times more massive on average. The two subsamples are statistically different in velocity 
width and Hi mass at 4.6cr- and 3.5(j-con£dence levels, respectively. Nearly every large 
spiral galaxy (i.e. IT 50 > 200km s“^, log Mhi > 9.5 Mq) detected by ALFA ZOA was also 
detected by 2MASS. The vast majority of smaller spirals and dwarf galaxies (i.e. W 50 < 
200 km s“\ log Mhi < 9.5 Mq) were not detected by 2MASS. 


6. Survey Performance 
6.1. Positional Accuracy 

The positional accuracy of ALFA ZOA should be well within the FWHM of the 
telescope (FWHM = 3.4') because the survey is Nyquist sampled. Figure fTOl shows the 
positional separation between ALFA ZOA Deep detections and their counterparts. The 
counterpart with the smallest separation is chosen for this plot when there are multiple 
possible counterparts for the same source. All of the nearest counterparts are 2MASS 
galaxies. The positional accuracy of 2MASS is 0.5" (Skrutskie et ah 2006), over two 
orders of magnitude finer than the FWHM of ALFA ZOA, meaning that the distribution 
of separations should be almost entirely due to the positional uncertainty of the ALFA 
ZOA Survey. There can be intrinsic offset between a galaxy’s neutral hydrogen and stellar 
structure, though an offset of a kiloparsec only subtends 3" at 4000 km s“^ (beyond which 
25 of 26 counterparts are located). 

Assuming that counterparts are real and the separations in Figure (TU] are due to 
uncertainty in ALFA ZOA positions, then their standard deviation is an estimate of the 
positional accuracy of the survey. The standard deviation of separations is cr = 0.47', and 
this value is adopted as the survey’s positional accuracy. The one (dashed circle) and 
two (solid circle) a boundaries are overplotted in Figure [TUI The distribution of positions 
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Fig. 10.— Separations (ALFA ZOA minus literature) in right ascension (Aa) and declination (Ai5) between 
ALFA ZOA Deep and NED counterparts, in arcminutes. The one (dashed circle) and two (solid circle) 
standard deviations of the separations are shown. 

shows an average offset in right ascension of < Aa >= 0.12' to the east. This offset is not 
considered to be statistically signihcant as it is fonnd to be at only a l.Sa-conhdence level. 
There is no indication that positions are systematically offset in declination. 
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6.2. Sensitivity 
6.2.1. Noise 

The noise level of the ALFA ZOA Deep data cube reached its expected rms value of 
1 mJy (at 9 km s“^ velocity resolution). Figure [TT] shows the cube’s rms as a function 
of heliocentric velocity. The rms was averaged over an inner quarter of the image plane, 
chosen in an area with a relatively low detection rate in order not to pollute the rms map 
with source flux. The mean noise is rms = 1 mJy for the high frequency sub-band (i.e. 
Vhei < 11,500 km s“^), but it increases to rms = 1.2 mJy for the low frequency sub-band 
(i.e. Vhei > 10,000 km s“^). The overlap between the two sub-bands can be seen between 
10,000 km s“^ and 11,500 km s“^. It is clear that the radio frequency interference (RFI) 
prevalent outside of the protected frequencies (i.e. 1400 - 1427 MHz) raises the average 
system temperature non-negligibly inside the low frequency sub-band. 
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Fig. 11.— Noise in mJy per heliocentric velocity averaged over the inner quarter of the image plane for 
both the high- and low-frequency sub-bands. Major RFI is labeled. 


Known RFI is labeled in Figure ITT] using information provided by Arecibo Observatory. 
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The GPS L3 satellite produces RFI at 1381 MHz that spans velocities from 8400 - 8800 km 
s“^. FAA radar from the airport at Punta Borinquen near Aguadilla and east of San Juan 
at Pico del Este produces RFI at 1350 MHz and 1330 MHz, covering velocities ranging from 
14,600 - 16,000 km s“^ and 19,700 - 20,700 km s“^, respectively. There is also RFI at 1339 
MHz, spanning velocities from 18,100 - 18,400 km s“^. There is no RFI source currently 
known to the Observatory at 1339 MHz. 


6.2.2. Hi Detection Limit 



Fig. 12.— Hi mass as a function of heliocentric velocity for both the ALFA ZOA Deep (open blue boxes) 
and Shallow (closed red triangles) Survey detections. Predicted Hi mass detection limit is shown for the 
Deep (dashed blue line) and Shallow (solid red line) Surveys. 


Figure [12] shows a plot of Hi mass as a function of velocity for both the Deep (open 
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boxes) and Shallow (closed triangles; rms = 5.4 mJy; Henning et al. 2010; McIntyre et al. 
in prep.) Surveys. The sensitivity and depth of the Deep Survey can clearly be seen. Gaps 
in detections at specihc redshifts are due to the RFI discussed above as well as large scale 
structure, which will be discussed in the next section. 

The detection limit of the survey is dehned as the farthest distance that a specihc H l 
mass can be detected above a signal-to-noise ratio, S/N = 6.5. A £t of the detection limit 
of each survey is shown in Figure [12] using criteria devised by Giovanelli et al. (2007) for a 
bivariate signal-to-noise ratio, 

^ rms (2 6 v IFso)^/"’ 

where Fhi is the integrated hux in Jy km s“^, IT 50 is the half peak velocity width in km 
s“^, rms is the noise in Jy, and Sv = 9 km s“^ is the velocity resolution of the survey. 
Though sources can be detected below the detection limit as seen in Figure Giovanelli 
et al. showed empirically that sources above S/N = 6.5 have a 95% reliability rate, and this 
estimate is adopted for this paper. Solving for Fhj and plugging into equation ([T|) gives an 
expression for the H l mass detection limit of the survey, 

Mhi = 2.36 X 10® S/N rms (2 Sv IF 5 o)^^^ (9) 

where r is the distance to a detection in Mpc, a value of H 50 = 200 km s“^ is chosen for the 
plot in Figure US] and the distance is converted from Mpc to heliocentric velocity using Hq 
= 70 km s“^. The typical noise of the survey is rms = 5.4 mJy, 1 mJy for the shallow, deep 
survey respectively. The superiority of the deep survey at detecting lower mass sources out 
to heliocentric velocities of 20,000 km s“^ can clearly be seen. 
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6.3. Completeness 

The completeness limit of the survey is the lowest integrated flux, Fnium, above 
which every galaxy can be detected. One technique for measuring completeness is by 
htting a power law with a slope of -3/2 to the histogram of flux. The completeness limit is 
reached where the histogram begins to deviate from the slope. However, this method only 
works if galaxies are homogeneous and isotropic, which is not the case for the small area 
surveyed in our First Results Catalog. An alternative method for determining the survey’s 
completeness is the Test for Completeness, Tc, a statistical test derived by Rauzy (2001) 
for a magnitude-redshift sample, independent of large scale structure. The method takes 
a sample of galaxies within a given volume and brighter than a given flux and compares 
the number of galaxies that are fainter and brighter than every galaxy in the sample. Tc is 
calculated for Hi as in Zwaan et al. (2004). 

Tc should follow a Gaussian distribution with an average value of 0 and unit variance 
for galaxy samples that are above the completeness limit, and move systematically to 
negative values for samples that are not complete. A plot of Tc as a function of integrated 
flux can be seen in Figure [T31 Tc begins to fall below -1 at Fniiim < 0.7 Jy km s h It 
hovers near -1 until Fnium < 0.4 Jy km s“^, where it systematically drops below -2. This 
indicates that the survey is not complete below 0.4 Jy km s“^ with about a 98% conhdence 
level. The completeness limit adopted for the ALFA ZOA Deep Survey is Fuiiim = 0.5 Jy 
km s“^, where Tc = -1.25 (i.e. the lowest value before the clear systematic drop in Tc). 

The bottom panel of Figure [13] shows an attempt to measure the completeness from 
where the flux histogram deviates from a -3/2 power law. The value for completeness 
using this method is Fuium = 1 Jy km s“^, signihcantly higher than the Tc method. As 
mentioned above, the power law technique should fail if the distribution of galaxies in 
the survey area is not homogeneous, which is the case for our small First Results area 
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Fig. 13.— Top Panel. Test of completeness for ALFA ZOA. Horizontal dashed lines indicate unit vari¬ 
ance. The vertical dashed line shows the adopted completeness limit. Bottom Panel. Log-log histogram of 
integrated flux. The -3/2 slope (diagonal dashed line) falls off from the fit somewhere between -0.1 and 0.1, 
indicating that the completeness is reached somewhere around 1 Jy km s“^ (vertical dashed line). 

dominated by large-scale structure. Further, a completeness limit of Fhi um = 0.5 Jy km s“^ 
makes sense compared to recent H l surveys, given differing rms sensitivities, as illustrated 
in Figure [TH 

For normalized comparison, the noise for ALFALFA (Saintonge 2007) and HIPASS 
(Zwaan et al. 2004) is scaled to a velocity resolution of 9 km s“^. Figure HH shows 
a remarkably well £t, linear relationship between Hi survey noise and integrated flux 
completeness limit, indicating that the Tc method is giving an expected value compared to 
how other surveys are performing. It should be noted that ALFALFA (Haynes et al. 2011) 











38 



Fig. 14.— Integrated flux completeness limit versus rms for major surveys. The open diamonds indicate 
surveys that used the Tc method and the closed circles are for surveys that fit a -3/2 power law for calculating 
completeness limit. 

published bivariate completeness as a function of integrated flux, Fhi, and galaxy velocity 
width, W 50 , and so the completeness limit for their average W^o is used for Figure [HI 

Figure [15] shows a log-log plot of integrated flux, Fhi, versus velocity width, IF 50 for 
ALFA ZOA. Similar to the bivariate signal-to-noise ratio used to estimate the detection 
limit of the survey, discussed in Section 16.2.21 the relationship between Fhi and W 50 can be 
used to estimate the bivariate completeness limit, FniiW^o) um, above which every galaxy 
can be detected. FniiW^o) um is estimated by scaling equation ([S]) so that Fhi um = 0.5 Jy 
km s“^ at the mean W^q. The mean velocity width is = 167 km s“^, so the bivariate 
completeness limit can be estimated as, 

log[FHiiW^o) Um] = 0.5 logiW^o um) — 1-4. (10) 

Both the bivariate detection limit (dashed line) and bivariate completeness limit (solid line) 
are shown in Figure [HI It should be noted that, while three ALFA ZOA Deep detections 
are located below the detection limit line in Figure [12] and none are below the line here, 
that is because an estimate of IF 50 = 200 km s“^ was assumed to make the plot in Figure 
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Fig. 15.— Log-log plot of integrated flux versus velocity width. The detection limit of the survey (dashed 
line) and the bivariate completeness limit (solid line) are shown. 


m 


6.4. Zone of Avoidance 

Figure [16] shows a plot of integrated flux, Fhj, as a function of foreground extinction, 
Ab- Extinction estimates were taken at the position of each ALFA ZOA detection using 
NED as described in Section l5dl It should be noted that estimates of extinction at |6| < 5° 
are particularly uncertain (Schlegel, Finkbeiner & Davis 1998). This plot does not show a 
relationship between Fhi and extinction, indicating that our survey does not have a Zone 
of Avoidance, as expected. 

Sources with higher Fhi are more likely to have a counterpart, and counterparts are 
less likely at higher Ab- This is expected, as the purpose of ALFA ZOA is to detect galaxies 
in areas with high extinction, where surveys at higher frequencies cannot. 55% of detections 
have a counterpart below Ab = 3.5 magnitudes while only 37% have counterparts above 
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Fig. 16.— Integrated flux as a function of extinction, Ab- Detections with a counterpart (open red circles) 
and with no counterpart (blue x’s) are shown. 

Ab = 3.5 mag. 

Because the number of detections is not very high, it is possible that the detections at 
high extinction happen to come from intrinsically low mass or low velocity width galaxies 
that are less sensitive to NIR observations. Figure ITTI shows a plot of extinction, Ab, as 
a function of Mhi and hFso, color coded for those with a counterpart (open red circles) 
and those without (blue x’s). There is a near uniform distribution of Mhi and IT 50 at all 
extinctions, indicating there is nothing intrinsic about the galaxies being detected that 
causes a lack of counterparts at high extinction. The even distribution of IF 50 and Mhi 
across extinction is another indicator of the efficacy of ALFA ZOA at penetrating the 
Zone of Avoidance. For instance, ALFA ZOA J0610+1709 has the highest Mhi, and 15th 
highest Fhi in the catalog, with a large velocity width W 50 = 215 km s“^, but it has no 
counterpart in the literature, even though its measured parameters indicate a massive spiral 
that should be easily detected by a deep NIR survey. The most likely reason why it has 
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Fig. 17.— Left panel. Extinction versus logMni, color coded for those with counterparts (open red circles) 
and those without (blue x’s). Right panel. Extinction versus velocity width, IF 50 , same color scheme as left 
panel. 

gone undiscovered is because it is located in an area of high extinction, Ab = 6.7 {Aj = 
1.3) or high stellar crowding, both of which affect low-latitude optical and NIR surveys but 
are of no concern to 21-cm surveys. 


7. Large Scale Structure 

While it is not quite a “pencil beam” survey the hrst results area is narrow enough 
that it does not fairly probe large scale structure (LSS) in angular dimensions. Therefore, 
examining a histogram of redshifts is a useful way to explore large scale structure in 
redshift space. Figure fT^ shows detections per velocity above the completeness limit (thick 
line) along with the expected number of detections per velocity (thin line). The expected 
number of detections is calculated by integrating the HIMF at the completeness limit 
using values from Martin et ah (2010). As we are expecting to detect every galaxy above 
the completeness limit, the difference between the two lines is a measure of the over- 
and underdensity of the large scale structure cut across by the survey. Redshifts where 
signihcant RFI occurs are labeled. RFI can limit the sensitivity of the survey over the 
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Fig. 18. — Left panel. Distribution of ALFA ZOA detections (thick solid line) that are above the complete¬ 
ness limit, as a function of redshift. The expected distribution (thin solid line) from integrating the HIMF 
(Martin et al. 2010) is shown along with a more robust expectation (dashed line) from taking into account 
RFI and rms = 1.2 mJy for the low sub-band. Velocities with significant RFl are labeled; the arrow with 
white polkadots indicates the redshift and bandwidth of the GPS L3 Satellite, black polkadots are for FAA 
radar. Right panel. The log difference between ALFA ZOA and the robust expected distribution. Observed 
large scale structure is shaded white. All lines are smoothed for visualization purposes using a cubic Bezier 
curve. 

volume affected by its bandwidth. A more robust estimate of the expected distribution of 
detections for ALFA ZOA is shown (dashed line) by taking into account the effect of RFI 
on survey volume as well as the increase in noise in the low sub-band (i.e rms = 1.2 mJy 
for Vhei > 11,500 km s“^). The lines in Figure fTSl are smoothed using a cubic Bezier curve 
for the sake of visualization. The smoothing is not used for analysis. 


7.1. Comparison to Predicted Large Scale Structure 

Erdogdu et al. (2006) created density reconstruction maps out to 16,000 km s“^ from 
the Two Micron Redshift Survey (2MRS; Huchra et al. 2012). The Zone of Avoidance for 
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2MRS is \b\ < 5° in the outer Galaxy, and so Erdogdu et al. were forced to predict the 
extent of large scale structure across the ZOA. Comparing ALFA ZOA detections with the 
Erdogdu density maps is a good check for the effectiveness of predicting LSS in the ZOA. 
Figure [19] shows a wedge plot of ALFA ZOA detections (open circles) alongside 2MRS 
detections (black dots) and indicates major overdensities from the density reconstructions 
of Erdogdu et ah inside the red ovals. ALFA ZOA confirms the continuation across the 
ZOA of much of the structure predicted by Erdogdu et ah, and contradicts some of the 
predictions as well. 

Referring to Figures [TS] and Hn] the prediction of the overdensity Cl and its extension 
through C5 from 2000 km s“^ to 4000 km s“^ is neither confirmed nor unconfirmed, though 
galaxies at the location of Cl and C5 are detected. The Orion overdensity at 6000 km s“^ 
is clearly conhrmed to extend into the ZOA, as is Abell 539 at 8000 km s“^. Abell 539 is 
most likely an even greater overdensity than detected here because of RFI from the GPS 
L3 satellite affecting velocities 8400 km s“^ - 8800 km s“^. Galaxies between Orion and 
Abell 539 can be seen around Vhei ~ 7000 km s“^, though whether there is a connection 
between Orion and Abell 539 in redshift space cannot be confirmed. CIZA J0603.8+2939 
is confirmed not to extend below 6 = 0°. A significant underdensity is confirmed between 
9000 and 12,000 km s“^, though this underdensity appears to be greater and extend farther 
in redshift than predicted. Large scale structure from C21 is confirmed as an overdensity 
starting at 13,500 km s“^ but its extent through 16,000 km s“^ cannot be confirmed 
because of FAA radar interference. The very strong RFI presence due to FAA radar from 
14,600 to 16,000 km s“^ most likely contributes to a lack of detections of C21 galaxies in 
that velocity range. ALFA ZOA detects an underdensity from 16,000 to 20,000 km s“^, at 
velocities beyond the range of 2MRS and the Erdogdu density reconstruction maps. There 
are no detections above the completeness limit in that velocity range though an average 
universe would contain 10 galaxies above the completeness limit in that region. This is a 
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Fig. 19.— 2MRS galaxies (black dots) and ALFA ZOA detections (blue open circles). Overdensities 
predicted from 2MRS (Erdogdu et al. 2006) are circled in red and the Erdogdu et al. naming convention is 
preserved. The field of view of the first results survey is shown (dotted line). In the plot, CIZA is shorthand 
for CIZA J0603.8+2939. 


statistically significant nnderdensity at abont a 99% confidence level. There is narrow RFI 
aronnd 18,500 km s“^ that affects sensitivity at that velocity and the FAA radar starting 
at 20,000 km s“^ rednces sensitivity to zero at the edge of the cube’s velocity cutoff. Figure 
[20] shows sky distribution plots of ALFA ZOA and 2MRS detections for the three volumes 
where ALFA ZOA detects an overdensity. While near infrared maintains a clear Zone 
of Avoidance, ALFA ZOA easily traces large scale structure through the lowest Galactic 
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latitudes. 


Orion Abell 539 C21 



Fig. 20.— Sky distribution plots showing 2MRS (black dots) and ALFA ZOA (open blue circles) detections 
for the three major overdensities detected by ALFA ZOA. The overdensities are labeled above each plot and 
the field of view of the first results survey is shown (dotted line). Orion covers the velocity range 5000 to 
6000 km s“^. Abell 549 ranges from 7500 to 8500 km s“^. C21 ranges from 12,000 to 15,000 km s“^. 


8. Conclusion 

This first data release of the ALFA ZOA Deep Survey demonstrates that the survey 
achieves its expected noise level of 1 mJy (at 9 km s“^ resolution) and it is complete above 
Fhi = 0.5 Jy km s“^. First results display the ability to detect galaxies out to nearly 
Vhei = 19,000 km s“^ and at extinctions that surveys at other wavelengths struggle to 
penetrate. ALFA ZOA Deep has completed 15 square degrees in the outer Galaxy and 
is continuing to take data, intending to cover about 300 square degrees in both the inner 
(30° < I < 75°; 6 < |2°|) and outer (175° < I < 207°; -2° < b < +1°) Galaxy, over the next 
several years. 
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As results from the full survey come in, asymmetry in large scale structure beyond 100 
Mpc, a possible explanation for the discrepancy in the mass density dipole compared to 
the CMB dipole, will be probed in regions of the inner and outer Galaxy. The ability to do 
so has been demonstrated by the survey’s first results. Also shown is a detection rate that 
indicates the full survey’s ability to detect over a thousand galaxies, which can be used to 
measure of the low mass end of the Hi Mass Function out to farther distances than have 
been measured previously by large, blind H l surveys. 
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Table 2. Catalog of possible counterparts for ALFA ZOA detections. 


ALFAZOA 

1 

n 

b 

n 

Ab 

(mag) 

Counterpart 

Sep. 

(') 

Av 

(km ) 

J0556 + 1450 

193.43 

-5.02 

1.2 

2MASX J05563099H-1450236 

0.3 

2 

J0556+1513 

193.14 

-4.74 

1.4 

2MASX J0556551bfl513206 

0.3 


J0556+1520 

192.95 

-4.84 

1.3 

2MASX J05561201+1520386 

0.5 






ZOAG G192.97-04.82 

1.3 






2MASX J05561519fl520186 

1.3 


J0558+1557 

192.73 

-3.98 

2.2 

2MASX J0558484&fl557068 

0.4 


J0558+1708 

191.61 

-3.53 

2.6 

2MASX J05581034H-1708233 

0.7 


J0559+1543 

193.02 

-3.95 

2.0 

2MASX J05593112+1542527 

0.8 






ZOAG G193.03-03.9 

0.9 


J0559+1636 

192.24 

-3.52 

2.9 

2MASX J05593052+1636117 

0.6 


J0601 + 1540 

193.25 

-3.66 

2.5 

2MASX J0601016C+1540102 

0.1 


J0600+1640 

192.32 

-3.21 

5.3 

2MASX J0600486C+1641193 

1.0 


J0600+1651 

192.20 

-3.11 

4.0 

2MASX J06005403fl649473 

1.4 






2MASX J06005104+1649293 

1.7 


J0600+1654 

192.12 

-3.13 

3.0 

2MASX J0600358bfl654227 

0.3 


J0601 + 1709 

191.97 

-2.89 

4.6 

2MASX J06011037+1709532 

0.6 


J0602 + 1438 

194.34 

-3.83 

3.2 

2MASX J06023859fl438273 

0.4 


J0602 + 1450 

194.13 

-3.77 

3.2 

2MASX J06022660+1450540 

0.2 

2 

J0603+1601 

193.19 

-3.06 

3.9 

2MASX J06030641+1601033 

0.3 


J0603+1607 

193.09 

-3.00 

4.2 

2MASX J06030454+1607413 

0.4 






Table 2—Continued 


ALFAZOA 

1 

n 

b 

n 

(mag) 

Counterpart 

Sep. Av 

(') (km ) 

J0605 + 1526 

194.03 

-2.75 

4.4 

2MASX J06055511+1526548 

0.8 

J0607+1506 

194.48 

-2.62 

3.6 

2MASX J06072476H-1505517 

1.9 

J0607+1543 

193.97 

-2.24 

7.0 

2MASX J06074014H-1544206 

0.7 

J0607+1607 

193.65 

-2.03 

5.4 

2MASX J06074379H-1608036 

1.1 

J0608+1527 

194.35 

-2.12 

5.3 

2MASX J06085110T1526540 

0.6 

J0611 + 1552 

194.25 

-1.44 

4.8 

2MASX J06110597+1551441 

0.5 

J0614+1641 

193.91 

-0.36 

7.4 

2MASX J06142201+1641578 

0.5 

J0616+1626 

194.39 

0.02 

5.9 

2MASX J06164466H-1627254 

0.5 

J0616+1641 

194.11 

0.02 

5.3 

2MASX J06161202+1641266 

0.7 

J0617+1655 

194.03 

0.36 

5.0 

2MASX J06171616H-1655575 

0.4 






